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There are many examples of wave equations in the physical sciences, char- 
acterized by oscillating solutions that propagate through space and time while, 
in lossless media, conserving energy. Examples include the scalar wave equa- 
tion (e.g. pressure waves in a gas). Maxwell's equations (electromagnetism) , 
Schrodingcr's equation (quantum mechanics), clastic vibrations, and so on. 
From an algebraic perspective, all of these share certain common features. They 
can all be written abstractly in a form 

dw 

- = Z.w + s (1) 

where w(x, t) is some vector-field wave function characterizing the solutions (e.g. 
a 6-component electric+magnetic field in electromagnetism), D is some linear 
operator (using the "hat" notation from quantum mechanics to denote linear 
operators), neglecting nonlinear effects, and s(x,i) is some source term. The 
key property of D for a wave equation is that it is anti-Hermitian, as opposed to a 
parabolic equation (e.g. a diffusion equation) where D is Hermitian and negative 
scmi-dcfinitc. From this anti- Hermitian property follow familiar features of wave 
equations, such as oscillating/propagating solutions and conservation of energy. 
In many cases, we will set s = and focus on the source-free behavior. 

In the following note, we first derive some general properties of cq. (1) from 
the characteristics of D, and then give examples of physical wave equations that 
can be written in this form and have these characteristics. 
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1 General properties of wave equations 

In lossless media, D turns out to be an anti-Hermitian operator under some 
inner product (w,w') between any two fields w{x,t) and w'(x, at a given 
time t. This means that (w,Z)w') = — (Dw.w') for any w, w': D flips sign 
when it moves from one side to the other of an inner product. This is proven 
below for several common wave equations. Formally, & = —D, where f is the 
adjoint: for any operator A, (w,4w') = (A^w,w') by definition of A'. 

The anti-Hermitian property of D immediately leads to many useful conse- 
quences, and in particular the features that make the equation 'Vave-like:" 

1.1 Harmonic modes and Hermitian eigenproblems 

First, we can write down an eigen- equation for the harmonic-mode solutions 
w(x, = W(x)e^"^', assuming D is linear and time-invariant. Substituting 
W(x)e~"^* into eq. (1) for the source-free s = case, we obtain: 

wW = ibW, (2) 

which is a Hermitian eigenproblem: if D is anti-Hermitian, then iD is Hermitian. 
It then follows that the eigenvalues u) are real and that solutions W can be 

chosen orthogonal. Notice that the real eigenvalues lo corresponds directly to our 
assumption that the medium has no dissipation (or gain) — if uj were complex, 
waves would exponentially decay (or grow), but instead they oscillate forever in 
time. For any reasonable physical problem, it also follows that the eigenmodes 
are a complete basis for all w, so they completely characterize all solutions.^ 

In a problem where D is purely real (pretty much every wave equation except 
for Schrodinger), the complex conjugate of an eigenfunction is also an eigenfunc- 
tion with the conjugated eigenvalue, so the eigenvalues come in ±iuj pairs — by 

-'Technically, we require D or its inverse to be compact, which is true as long as we have 
a reasonable decaying Green's function This is almost always true for physical wave 

problems, but functional analysts love to come up with pathological exceptions to this rule. 



2 



adding complex conjugates, one can therefore get the purely real solutions one 
expects in a real equation. 

1.2 Plane wave solutions 

Second, the most familar feature of wave equations is the existence of not just 
oscillation in time (real uj). hut oscillation in space as well. In particular, with 
a wave equation one immediately thinks of sinusoidal planewave solutions ~ 
gj(k x-wt) gQj^g iijfive vector k. These solutions arise in any equation of 
the anti-Hermitian form (1) that additionally has translational symmetry (the 
medium is homogeneous). 

Intuitively, translational symmetry means that D is the same at different 
points in space. Formally, we can make this precise by defining a translation 
operator Td that takes a function w(x,t) and translates it in space by a dis- 
placement d: 

fdw(x,t) = w(x-d,t). 

When we say that D has translational symmetry, we mean that D is the same if 
we first translate by some d, then operate D, then translate back: D = T^^DT^, 
or equivalently D commutes with Tj: 

Dfd = fat) 

for all displacements d. When one has commuting operators, however, one can 

choose simultaneous eigenvectors of both operators [2] . That means that the 
eigenvectors £'(x) of iD (and D) can be written in the form of eigenfunctions 
of Td, which are exponential functions e^^'^ for some k: 

w(x, t) = W{x)e-"^* = Woe'^'"-^-'"*) (3) 

for some constant vector Wq (depending on k and w) determined by D. The 
wave vector k must be real if we require our states to be bounded for all x (not 
exponentially growing in any direction).^ For each k, there will be a discrete set 
of eigenfrequencies cj„(k), called the dispersion relation or the band structure of 
the medium. 

Let us also mention two generalizations, both of which follow from the 
broader viewpoint of group representation theory [3, 4]. First, the existence 
of planewave solutions can be thought of as a consequence of group theory. 
The symmetry operators that commute with D form the symmetry group (or 
space group) of the problem (where the group operation is simply composi- 
tion), and it can be shown that the eigenfunctions of D can be chosen to trans- 
form as irreducible representations of the symmetry group. For the translation 

■^Complex k solutions are called evanescent waves, and can appear if D is only transla- 
tionally invariant over a finite region (e.g. if the medium is piecewise-constant). Even the 
real k solutions are a bit weird because they are not normalizable; ||e*'''^|p is not finite! By 
admitting such solutions, we are technically employing a "rigged" Hilbert space, which requires 
a bit of care but is not a major problem. Alternatively, we can put the whole system in a 
finite L x L x L box with periodic boundary conditions, in which case the components of k 
are restricted to discrete multiples of 27r/L, and take the limit L — » oo at the end of the day. 
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group {Td I d G K'^}, the irreducible representations are the exponential func- 
tions {e"*'''''}, but more complicated and interesting representations arise when 
one includes rotations and other symmetries. Second, in order to get wave so- 
lutions, one need not require D to commute with for all d. It is sufficient to 
require commutation for d = R = riiRi + n2R2 + n^Rs on a discrete periodic 
lattice R with primitive lattice vectors R^ (and any ne gZ): D is periodic, with 
discrete translational symmetry. In this case, one obtains the Bloch-Floquet 
theorem (most famous in solid-state physics): the eigenfunctions W can be cho- 
sen in the form of a plane wave multiplied by a periodic Block envelope. More 
explicitly, one has Block wave solutions: 

w(x,t) = Wk,„(x)e'['"^-'"" (4) 

where Wk,n(x) is a periodic function (invariant under translation by any R in 
the lattice) satisfying the Hermitian eigenproblem: 

The planewave case (3) [continuous translational symmetry] is simply the special 
case of the Bloch wave (4) [discrete translational symmetry] in the limit where 
the lattice vectors become infinitesimal (|R^| ^0). 



1.3 Time evolution and conservation of energy 

Third, we obtain conservation of energy in the absence of sources (s = 0), where 
we define the "energy" of a field w as its norm ||w|p = (w,w). The proof that 
II w IP is conserved for s = (sources would add or remove energy) is elementary, 
given an anti-Hermitian D: 

^|I.^|I2 Q ^ ^ ^ ^ 

— — — = — (w,w) = (w, w)-|-(w, w) = {Dw,w)+{w, Dw) = (Dw, w) — (Dw, w) = 0. 

This works even if D is time-dependent, which may seem surprising: if you take 
a wave equation and "shake it" by varying D rapidly in time, you might think 
you could add energy to the system. But no: a time-varying D can alter the 
frequency of the solution (which is not conserved in a time- varying problem), but 
not the energy. However, this is not the whole story, because D is not the only 
possible source of time-dependent behavior: the definition of our inner product 
(w, w') can depend on t as well. In fact, we will see that this is often possible for 
physical systems such as Maxwell's equations or even the scalar wave equation. 
In particular, our inner product is often of the form (w, w') = (w, Pw')o, where 
(•, )o denotes an inner product independent of time and P is some positive- 
definite Hermitian operator depending on the wave medium, which may depend 
on time. In this case, when taking the time derivative of ||w|p, we also get a 

term (w, ^w)o, which is not in general zero. So, a time-varying medium can 
break conservation of energy if the time variation ckanges tke norm. 
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If D is time-independent, we can easily write down the explicit solution of 
the initial-value problem. In this case eq. (1) for s = is solved formally by the 
operator exponential: 

w(x,i) = e'^*w(x,0) = C/tw(x,0) 

for an initial condition w(x, 0) and a time-evolution operator Ut = e^*. Because 
D is anti-Hermitian, it flips sign when it switches sides in an inner product, and 
hence Ut changes from Ut = e^* to (/} = e~^* = If^^ = tj-t- This means that 
tft is a unitary operator, and hence 

2 , „ , _ , 2 

[/jw = [UtW,Utw) = {Uf C/tW,w) = ||w|| . 

Thus, ||w|p does not change as the fleld evolves in time: energy is conserved! 

1.4 OfF-diagonal block form and reduced- rank eigenprob- 
lems 

There is another variant on the eigenequation (2) tjW — iUW: if we operate 
iD again on both sides, we get w^W = — i)^W, where —D^ is automatically 
Hermitian and positive semi-definite. Why would we want to do this? The main 
reason is that D often (not always) has a special block form: 

(5, 

for some operators D\ and D2. The anti-Hermitian property D* = —D implies 

that D2 = —Di and £)* = — -D2 under an appropriate inner product (see below). 
In this case, is block diagonal and negative-semidefinite: 



D2D1 \ / -DlDi 
D1D2 / V -D*2D2 



This means that the w^W = — D^W problem breaks into two lower-dimensional 
eigenproblems with operators ~D2Di — D^Di and —D1D2 = D2^2- In par- 
ticular, let us break our solution w into two pieces: 

Wi 

where Dk operates on Wfe, and suppose that our inner product also divides 
additively between these two pieces: 

(W,w') = (Wi,wi)i + (W2,W2)2 

in terms of some lower-dimensional inner products (•,•)! and (•, •)2- In this case, 
it immediately follows that —D2D1 = D\Di is Hermitian positive semi-definite 
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under (-j-)! and —D1D2 = D2D2 is Herniitian positive semi-definite under 
(•, •)2. We therefore have obtained two smaller Hermitian positive semi-definite 
eigenproblems 

-D2D1W1 = w^Wi, (6) 

-D1D2W2 = w^Wa, (7) 

again with real oj solutions and orthogonality relations on the W^..'^ 

Moreover, each of these has to give all of the eigenfrequencies co. Every 
eigenfunction of D must obviously be an eigenfunction of D^, and the converse 
is also true: given an eigenvector Wi of cq. (6) with eigenvalue w^, to get an 
eigenvector W of Z) with eigenvalue ui we just set 

W2 = -DiWi. (8) 

CO 

This formula looks a bit suspicious in the case where u) = 0: the static (non- 
oscillatory) solutions. For these (usually less-interesting) static solutions, the 
Wi and W2 eigenproblems decouple from one another and we can just set 
W2 = in W. Similarly, if we solve for an eigenvector W2 of cq. (7), we can 
construct W via Wi = iD2SN'2/uj for w 7^ or set Wi = for l^,' = 0. 

We lost the sign of uj by squaring it (which is precisely why we can solve 
an eigenproblem of "half" the size), but this doesn't matter: oj ~ ±^/uj^ both 
yield eigensolutions W by (8). Therefore, in problems of the block form (5), 
the eigenvalues lo always come in positive /negative pairs. In the common case 
where Di and D2 are purely real, the eigenvectors Wi (or W2) can also be 
chosen real, we can therefore obtain real solutions w by adding the +u! and —uj 
eigenfunctions (which are complex conjugates). 



1.5 Harmonic sources and reciprocity 

The most important kind of source s is a harmonic source 

s(x,i) = S(x)e-''"*. 

In this case, we are looking for the steady-state response w = W(x)(3^*'^*. 
Substituting these into eq. (1) problem of solving for W(x) is now in the form 
of an ordinary linear equation, rather than an eigenproblem: 

{ib - W = 6W = -iS. 

Notice that the operator O = iD — u on the left-hand side is Hermitian. Suppose 
that we solve the equation twice, with sources S'^^' and S^^^ to get solutions W^-'^^ 

•^The fact that they are positive semi-definite (and often positive-definite, since it is typi- 
cally possible to exclude the a; = solutions) is especially advantageous for numerical methods. 
Some of the best iterative eigensolver methods are restricted to positive-definite Hermitian 
problems, for example [5]. 
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and W*^^^. Then, the Hermitian property means that we obtain the following 
identity for the inner product: 

(ww,s(2)) = (w(i),i6w(2)) = (-z6w(i),w(2)) = -(s(i),w(2)). 

This is almost the same as a very well known property of wave equations, known 
as reciprocity. The reason it is only almost the same is that reciprocity relations 
normally use an unconjugated "inner product," assuming O is not only Hermitian 
but real (or complex-symmetric). This gets rid of the minus sign on the right- 
hand side, for one thing. It also only requires iD to be complex-symmetric 
(Hermitian under an unconjugated inner product) rather than Hermitian, which 
allows reciprocity to apply even to systems with dissipation. (It also simplifies 
some technical difficulties regarding the boundary conditions at infinity.) See 
e.g. Ref. [6].^ 

What if O is not real or complex-symmetric, e.g. in the common case where 
D is real-antisymmetric? Can we still have reciprocity with an unconjugated 

inner product? In this case, we can often instead express reciprocity using the 
block form (5), since in that case the operator is real-symmetric if D is real, as 
described below. 

Just as for the eigenproblem, it is common when we have the block form 
(5) of D to break the problem into a smaller linear equation for Wi or W2, 
similarly subdividing S = (Si;S2). For example, the equation for Wi is 

{-D2D1 - a;2)Wi = -iwSi + D2S2. (9) 

For example, in the case of the scalar wave equation below, this equation be- 
comes a inhomogeneous scalar Helmholtz equation. Again, notice that the op- 
erator —D2D1 —LiP' on the left-hand side is Hermitian. In fact, it is often purely 
real-symmetric (or complex-symmetric in a system with dissipation) , which al- 
lows us to derive recprocity using an unconjugated inner product as described 
above. 



2 The scalar wave equation 

There are many formulations of waves and wave equations in the physical sci- 
ences, but the prototypical example is the (source- free) scalar wave equation: 



where u(x, t) is the scalar wave amplitude and c = v ab is the phase velocity of 
the wave for some parameters a(x) and &(x) of the (possibly inhomogeneous) 
medium."* This can be written in the form of eq. (1) by splitting into two first- 
order equations, in terms of u and a new vector variable v satisfying aVu = v 
and bV -v = ii: 

dw d / u \ ( 6V- 



dt dt \ y J \ aV / V V 



u 



Dw 



''More generally, a may be a 3 x 3 tensor in an anisotropic medium. 
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for a 4 X 4 linear operator D and a 4-component vector w = (u; v), in 3 spatial 
dimensions. 

Next, we need to show that D is anti-Hermitian, for the case of lossless media 

where a and b are real and positive. To do this, we must first define an inner 
product (w,w') by the integral over all space: 

(w,w') = j [u* {b-^u')+w* {a-W)] d^x, 

where * denotes the complex conjugate (allowing complex w for generality). 

Now, in this inner product, it can easily be verified via integration by parts 
that: 

(u, bn') = J [u*V • v' + V* • Vu'] d^x, = ■■■ = -(£)u,u'), 

which by definition means that D is anti-Hermitian. All the other properties — 
conservation of energy, real eigenvalues, etcetera — ^then follow. 

Also, note that the —D eigenproblem in this case gives us more-convenient 
smaller eigenproblems, as noted earlier: — 6V- (aVw) = u^u, and — aV(6V- v) = 
u^v. (These may not look Hermitian, but remember our inner product.) And 
in the harmonic-source case, we get an operator —D2D1 — = — &V • aV — 
when solving for u via (9) — for the common case a = b = 1, this is the scalar 
Helmholtz operator — (V^ -|- w^). 



3 Maxwell's equations 

Maxwell's equations, in terms of the electric field E, magnetic field H, dielectric 
permittivity e and magnetic permeability /i, are (neglecting material dispersion 
or nonlinear effects): 

5(mH) 



V X E = 



dt 



where J is a current source; J = in the source- free case. The other two 
Maxwell's equations V • eE = and V • = express constraints on the fields 
(in the absence of free charge) that are preserved by the above two dynamical 
equations, and can thus be ignored for the purposes of this analysis (they are 
just constraints on the initial conditions). For simplicity, we restrict ourselves 
to the case where e and ^ do not depend on time: the materials may vary with 
position X. but they are not moving or changing."'' e and /i may also be 3x3 
tensors rather than scalars, in an anisotropic medium, but we assume that they 
are Hermitian positive-definite in order to be a lossless and transparent medium. 
In this case, the equations can be written in the form of (1) via: 



dw _ d f E \ _ f iVx \ f E \ f -J/s 



at dt\ B. J V -^vx J \ n J ' \ 



One can get very interesting physics by including the possibility of time- varying materials! 
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where w here is a 6-coniponent vector field. 

To show that this D is anti-Hermitian, we must again define an inner product 
by an integral over space at a fixed time: 

(w, w') = ^ /[E* • (£E') + H* • (mH')] d^x, 

which is precisely the classical electromagnetic energy in the E and H fields (for 
non-dispersive materials) [7]! 

Given this inner product, the rest is easy, given a single vector identity: for 
any vector fields F and G, V • (F x G) = G • (V x F) - F • (V x G). This vector 
identity allows us to integrate Vx by parts easily: / F • (V x G) = / G • (V x F) 
plus a surface term (from the divergence theorem) that goes to zero assuming 
||F|| and ||G|| are < oo. 

Now we can just plug in (w, DW) and integrate by parts: 

(w, ^w') = ^ y"[E* • (V X H') - H* • (V X E')] d^x 

= ^ y"[(V X E)* • H' - (V X H)* • E'] d^x = (-£)w, w'). 

Thus, D is Hermitian and conservation of energy, real uj, orthogonality, etcetera 
follow. 

Again, D is in the block form (5), so the eigenproblem simplifies into 
two separate Hermitian positive semi-definite eigenproblems for E and H: 

^Vx (^VxE) = a;2E, (11) 

^Vx (ivxH) = w2jj_ (12) 

These convenient formulations are a more common way to write the electro- 
magnetic eigenproblem [8] than (2). Again, they may not look very Hermitian 
because of the 1/e and 1/ ji terms multiplying on the left, but don't forget the 
e and /i factors in our inner product, which makes these operators Hermitian.^ 
Therefore, for example, two different eigensolutions Ei and E2 are orthogonal 
under the inner product (Ei,E2)e = JE* • {£^2) = 0. Bloch's theorem (4) for 
this Hermitian eigenproblem in the case of periodic materials leads to photonic 
crystals [8]. 

Recall that something funny happens at w = 0, where E and H decouple. In 
this case, there are infinitely many a; = static-charge solutions with V-eE 7^ 



"Alternatively, if we defined inner products without the e and /i, we could write it as 
a generalized Hermitian eigenproblem. For example, if we used the inner product (E, E') = 
f E* 'E', we would write the generalized Hermitian eigenproblem V X V X E) = o^^eE. This 

gives the same result: for generalized Hermitian eigenproblems Au = XBu, the orthogonality 
relation is (ui,Bu2) = for distinct eigenvectors ui and U2. This formulation is often more 
convenient than having s or n factors "hidden" inside the inner product. 
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and/or V • /iH ^ 0. For w ^ 0, V • eE = and V • /iH = follow automatically 
from the equations — iweE = V x H and iw/uH = V x E. (These divergence 
equations are just Gauss' laws for the free electric and magnetic charge densities, 

respectively.) 

In the case with a current source J ^ 0, we can solve for E via eq. (9), which 
here becomes: 

( - V X -V X -w^ ) E = iwJ/e, 
\e \x ) 

or (multiplying both sides by e) 

X X -b?e^ E = iwJ. 

The operator on the left side is Hermitian in lossless media, or complex-symmetric 
in most dissipativc materials where e and /i arc complex scalars. This gives rise 
to the well-known Rayleigh-Carson and Lorentz reciprocity relations [9, 6]. 



4 The one-way scalar wave equation 

As a break from the complexity of Maxwell's equations, let's look at the source- 
free one-way scalar wave equation: 

du du 
dx dt ' 

where c(x) > is the phase velocity. This is a one-way wave equation because, 
for c = constant, the solutions are functions u{x,t) = f(x — ct) traveling in the 
+x direction only with speed c. 

This equation is already in our form (1), with D = — c^, which is obviously 
anti-Hermitian under the inner product 

, p u*u' 
{u,u ) — / ax, 

J-oo C 

via elementary integration by parts. Hence energy is conserved, we have real 
eigenvalues, orthogonality, and all of that good stuff. We don't have the block 
form (5), which is not surprising: we can hardly split this trivial equation into 
two simpler ones. 



5 The Schrodinger equation 

Perhaps the most famous equation in form (1) is the Schrodinger equation of 
quantum mechanics: 

d'ip if fi^ ^2 tA , 



dt h \ 2m 
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for the wave function ')Jj{x,t) of a particle with mass m in a (real) potential 
y(x). Famous, not so much because Schrodinger is more well known than, 
say. Maxwell's equations, but rather because the Schrodinger equation (unlike 
Maxwell) is commonly taught almost precisely in the abstract form (1) [2, 10], 
where the operator D is given by 

for the Hamiltonian operator H. This is the equation where many students 
delve into abstract Hermitian operators and Hilbert spaces for the first time. 

The fact that H is Hermitian under the inner product {ij}, ij)') = J il>*'4>', and 
hence D is anti-Hermitian, is well known (and easy to show via integration by 
parts) . From this follow the familar properties of the Schrodinger equation. The 
eigenvalues w are interpreted as the energies E = hw, which are real (like any 
quantum observable), and the eigenstates ip are orthogonal. Conservation of 
"energy" here means that / l^/ip is a constant over time, and this is interpreted 
in quantum mechanics as the conservation of probability. 

Again, we don't have a block form (5) for the Schrodinger problem, nor do 
we want it. In this case, the operator D is not real, and we do distinguish 
positive and negative lj (which give very different energies E\). We want to 
stick with the full eigenproblem Dtjj = loiJj, or equivalently HiIj = Eip, thank 
you very much. Bloch's theorem (4) for this Hermitian eigenproblem in the case 
of a periodic potential V leads to the field of solid-state physics for crystalline 
materials [11, 12]. 

6 Elastic vibrations in linear solids 

One of the more complicated wave equations is that of elastic (acoustic) waves 
in solid media, for which there are three kinds of vibrating waves: longitudinal 
(pressure/compression) waves and two orthogonal transverse (shear) waves. All 
of these solutions are characterized by the displacement vector u(x,t), which 
describes the displacement u of the point x in the solid at time t. However, 
it will turn out that u itself is not the natural dynamic variable with which 
to characterize the wave equation, because it does not yield our preferred anti- 
Hermitian form in which energy-conservation is obvious. To begin with, we 
will consider the common simplified case of isotropic solids, and generalize to 
the anisotropic case in Sec. 6.4. We will neglect material dispersion (a valid 
approximation in a narrow enough bandwidth for nearly lossless materials) and 
consider only the linear regime of sufficiently small displacements. 

6.1 Isotropic materials 

In the linear regime, an isotropic elastic medium is characterized by its density p 
and the Lame elastic constants /j, and A: fi is the shear modulus and X = K—^fj,, 
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where K is the bulk modulus. All three of these quantities are functions of x for 
an inhomogeneous medium. The displacement u then satisfies the Lame-Navier 
equation [13, 14]: 

= V(AV • u) + V • [m (Vu + (Vu)^)] + f , (13) 

where f (x, i) is a source term (an external force density). This notation requires 
a bit of explanation. By Vu, we mean the rank-2 tensor (3x3 matrix) with 
(m, n)th entry dum/dxn (that is, its rows are the gradients of each component 
of u), and (Vu)-^ is the transpose of this tensor. By the divergence V- of such 
a tensor, we mean the vector formed by taking the divergence of each column 
of the tensor (exactly like the usual rule for taking the dot product of a vector 
with a 3 X 3 matrix). 

Clearly, we must break eq. (13) into two first-order equations in order to 
cast it into our abstract form (1), but what variables should we choose? Here, 
we can be guided by the fact that we eventually want to obtain conservation of 
energy, so we can look at what determines the physical energy of a vibrating 
wave. The kinetic energy is obviously the integral of |p|vp, where v is the 
velocity li, so v should be one of our variables. To get the potential energy, is 
convenient to first define the (symmetric) strain tensor s [14]: 

e=l{Vu+iVnf), (14) 

which should look familiar because 2e appears in eq. (13). In terms of £, the 
potential energy density is then ij,tv{e^s) + A| tr £p/2, where f is the conjugate- 
transpose (adjoint) [14]. Note also that tre = V • u, by definition of e. 

Therefore, we should define our wave field w as the 9-component w = (v; e), 
with an inner product 

(w, w') = [pv* ■ v' + 2^tr(£t£') + A(tr£)*(tr£')] d^x, 

so that jjwjj^ is the physical energy. Now, in terms of v and £, our equations of 
motion must be: 



dt dt\ e \ b, j I £ j ^ I 



£)w + s, (15) 



with 



£ = Div = ^ (Vv + (Vv)^) , (16) 

v = £)2£= i[V(Atr£)+2V-(M£)]. (17) 

Now, we can check that D is anti-Hermitian (assuming lossless materials, i.e. 
real A and fj). For compactness, we will use the Einstein index notation that 
repeated indices are summed, e.g. a • b = anbn with an implicit sum over n 
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since it is repeated, and let dk denote Then, plugging in the definitions 

above and integrating by parts: 



(w,Z)w') = - / pw* ■ {D2e') + 2/itr(£t£)iv') + A(tre)*(tr£)iv') 



= (I)w,w'). 

Note that we used the symmetry of e to write 2tv{e'^ Diw') = Smn^mv'n + 

We are done! After a bit of effort to define everything properly and prove 
the anti-Hermitian property of D, we can again immediately quote all of the 
useful general properties of wave equations: oscillating solutions, planewaves in 
homogeneous media, conservation of energy, and so on! 



6.2 Scalar pressure waves in a fluid or gcis 

For a fluid or gas, the shear modulus ji is zero — there are no transverse waves — 
and A is just the bulk modulus. In this case, we can reduce the problem to 
a scalar wave equation in terms of the pressure P = — AV • u = — tre. We 
obtain precisely the scalar wave equation (10) if we set a = 1/p, 6 = A, m = P, 
and V = — li (the velocity, with a sign flip by the conventions defined in our 
scalar- wave section). The wave speed is c = i/A/p. Our "energy" in the scalar 
wave equation is again interpreted as (twice) the physical energy: the integral 
of (twice) the potential energy P^/A and (twice) the kinetic energy />|v|^. 



6.3 Eigenequations and constraints 

Again, we have the block form (5) of D, so again we can write an eigenproblem 
for (J^ by solving for harmonic v or e individually. Given a harmonic u = 

U(x)e~*'^*, the equation for U just comes from the eigen-equation for v = 
— iwUe"*'^*, which is just the same as plugging a harmonic u into our original 
equation (13): 

w^U = - ^ {V(AV • U) + V • (VU + (VU)"^)] } , 

which again is a Hermitian positive semi-deflnite eigenproblem thanks to the 
factor of p in our inner product (U,U') = j pU* • U'. (Alternatively, one 
obtains a generalized Hermitian eigenproblem for an inner product without p.) 
Bloch's theorem (4) for this Hermitian eigenproblem in the case of periodic 
materials leads to the study of phononic crystals [15]. 
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We could also write the eigenequation in terms of e, of course. Sometimes, e 
is more convenient because of how boundary conditions are expressed in elastic 
problems [14]. However, in the case of e we need to enforce the constraint 

V X (V X e)^ = 0, (18) 

which follows from the definition (14) of e in terms of u. Note that the curl of a 
tensor, here, is defined as the curl of each column of the tensor, so V x (V x e)^ 
means that we take the curl of every row and column of e. This is zero because 
e is defined as the sum of a tensor whose rows are gradients and its transpose, 
and the curl of a gradient is zero. The same constraint also follows from eq. (16) 
for aj 7^ 0, which is part of the eigenequation, but e is decoupled from eq. (16) 
for w = and in that case we need to impose (18) explicitly [14]. 

6.4 Anisotropic linear materials 

The above equations assumed that the solid is isotropic, e.g. that it responds in 

the same way to compression along any axis. The more general equation involves 
a rank-4 elastic-modulus tensor C (a "four-dimensional matrix") with compo- 
nents Ciji^f, [16]. This tensor must satisfy the symmetries Cj^fe^ = Cjiki = Ckuj 
(implied by time-reversal symmetry [17] and also by more basic considera- 
tions [18], and furthermore C must be real in the absence of dissipation), which 
means that it is characterized by 21 parameters [16]. 

The definition (14) of the strain tensor e is the same, and we still write w = 
(v, e) in terms of the same dynamical variables v and e and the same off-diagonal 
form (15) of the wave equation. The relationship (16) between e and v is 
unchanged, so Di is unchanged. However, D2 and the stress-strain relationship 
which gives equation (17) for v (via Newton's law) is now generalized to: 

v = £)2£= •(€:£), (19) 

where by C : £ we mean the tensor contraction of the last two indices of C 
with those of e. In Einstein index notation (implicit sums of repeated indices) 
and denoting by di as above, eq. (19) is equivalent to Vi = djCjikiSu = 

Now, the key to showing that this wave equation is still anti-Hermitian is to 
choose the correct inner product: 

(w,w') = ^ j [pv*-w + e*:C:e']<f^, 

where by : C : e' we mean the contraction e*^Cijke£'ki = {^'ik^kUj£ji)* , which 
still has the correct symmetry (thanks to the symmetry of C) to make this a 
proper inner product. This is chosen to give the correct expression for the energy 
(w, w) = I /[p|v|^ + e* : C : e] [16]. It is now easily verified that = -D as 
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desired, and hence energy is conserved: 



(w, £)w') 



1 



pv* • {D2e') + e* :C: {Di^') 



dkv'e + dtv'u 



= \ j [vtdjCjikie'f,e + e*^C,jkidkv't,] (fix 

= -^y [s'keCkijidjV* + v'^dkCijkie*i\ cfi: 
= {-Dw,w'). 



d'^x 



7 The scalar wave equation in space 

Now, let's consider something a little different. We'll start with the scalar wave 
equation from above, aVu = v and 6V ■ v = ii, but instead of solving for u 
and v as a function of time t, we'll solve the equations as a function of space 
for harmonic modes u{x,t) = [/(x)e~"^* and v(x,t) = V(x)e^"^''' with uj ^ 0. 
In particular, we'll pick one direction, z, and look at propagation along the z 
direction. Plugging these harmonic modes into the scalar wave equations, and 
putting all of the d/dz derivatives on the left-hand side, we obtain: 



dw 



- _^ 
dz 

'" dx 



u 







-5- + Vt 



u 



Dw, 



y denotes the "transverse" (xy) del operator, and we have 



where 

eliminated the transverse components Vt of V via —iLuVt = aVfU. This is 
again of the form of our abstract equation (1), with z replacing t\ But is this D 
anti-Hermitian? 

Let's define our "inner product" via an integral over the xy plane: 



(w,w')= J [U*V^ + U'V:] d^x. 



(20) 



This is actually not a true inner product by the strict definition, because it is 
not positive definite: it is possible for ||w|p = (w,w) to be non-positive for 
w ^ 0. We'll have to live with the consequences of that later. But first, let's 
verify that (w,£>w') = — (£)w, w') for real a and 6, via integration by parts: 



(w,Z)w') 



u*[-^u'] + u*Vf-VtU'+(--v^]v: 

b I iw \ a 



d'x 



-{Dw,w'). 



d'x 
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So D is "anti-Hcrmitian," but in a "fake" inner product. To figure out the 
consequences of this, we have to go back to our original abstract derivations 
and look carefully to see where (if at all) we relied on the positive-definite 
property of inner products. 

Looking back at our proof, we see that we didn't rely on positive-definiteness 
at all in our derivation of conservation of energy; we just took the derivative of 
(w, w). So, "energy" is conserved in space. What docs that mean? Equation (20) 
can be interpreted as an integral of a time-average energy flux through the xy 
plane at z. Conservation of energy means that, in the absence of sources, with 
a steady-state (harmonic) solution, energy cannot be building up or decaying 
at any z. 

What about the Hermitian eigenproblem? What is the eigenproblem? We 

already have time-harmonic solutions. We can set irp an eigenproblem in z only 
when the problem (i.e., a and b) is z-invariant (analogous to the time-invariance 
required for time-harmonic modes). ^ In a 2;-invariant problem, the z direction 
is separable and we can look for solutions of the form •w{x,y,z) = W(a;,y)e'^^, 
where /3 is a propagation constant. Such solutions satisfy 

WW = (3W, 

which is a Hermitian eigenproblem but in our non-positive norm. Now. if one 
goes back to the derivation of real eigenvalues for Hermitian eigenproblems, 
something goes wrong: the eigenvalue /3 is only real if (W, W) ^ 0. 

This, however is a useful and important result, once we recall that (W,W) 
is the "energy flux" in the z direction. Modes with real /? are called propagating 
modes, and they have non-zero flux: propagating modes transport energy since 
(W,W) ^ for a typical real-/? mode.^ Modes with complex (most often, 
purely imaginary) /3 are called evanescent modes, and they do not transport 
energy since (W, W) = for them. 

What about orthogonality? Again, going back to the derivation of orthogo- 
nality for Hermitian eigenproblems, and not assuming the eigenvalues are real 
(since they aren't for evanescent modes). The orthogonality condition follows 
from the equation [/3* — (3'] ■ (W, W) = for eigensolutions W and W with 
eigenvalues /3 and /?'. That is, modes are orthogonal if their eigenvalues are 
not complex conjugates. Moreover, since iD is purely real, the eigenvalues must 
come in complex-conjugate pairs. This very useful, because it means if we want 
a fleld "parallel" to some eigenfield W and orthogonal to all the other eigenflelds, 
we just use the eigenfleld with the conjugate eigenvalue to W. And since D is 
real, that is just W*! So, in short, we have: 

(W*,W') =0 

for distinct eigenvectors W and W: the eigenvectors are orthogonal under the 
unconjugated inner product." 

'^Actually, we only need it to be periodic in z, in which case we can look for Bloch modes, 
is possible to have standing-wave modes with real /3 and (W, W) = 0, but these mainly 
occur at /9 = (where symmetry causes the flux to be zero). 
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Why have we gone through all this analysis of what seems like a rather 
odd problem? Because this is a simpler version of something that turns out to 
be extremely useful in more complicated wave equations. In particular, the z- 
propagation version of the time-harmonic Maxwell equations plays a central role 
in understanding waveguides and "coupled-wave theory" in electromagnetism 
[19]. 
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